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Abstract 

For a two-component Fermi gas in the unitarity limit (ie, with infinite scattering 
length), there is a well-known virial theorem, first shown by J. E. Thomas et al, 
Phys. Rev. Lett. 95, 120402 (2005). A few people rederived this result, and extended 
it to few-body systems, but their results are all restricted to the unitarity limit. Here 
I show that there is a generalized virial theorem for FINITE scattering lengths. I also 
generalize an exact result concerning the pressure, first shown in cond-mat/0508320, 
to the case of imbalanced populations. 
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Two-component ultracold atomic Fermi gases with large scattering lengths 
have been realized recently, and have become a focus of numerous research 
activities. In this paper we study two exact properties of such a system. 

We will consider the zero-range interaction model only, in which the scattering 
length a between the f and [ spin states is the only parameter for the inter- 
action. (Such a model is justified by typical experimental setups, in which the 
interatomic distance, the thermal de Broglie wavelength, and a are all large 
compared to the Van de Waals range of the interaction.) 



1 Generalized Virial Theorem 



If the system is confined by a harmonic trap, and is in the unitarity limit 
(a — » oo), the total energy E is related to the external potential energy by 



E = 2E 
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The result ([I]) was first shown by J. E. Thomas et al using the local density 
approximation pQ. A few people rederived this result, and extended it to few- 
body systems, but their results are all restricted to the unitarity limit and to 
a harmonic confinement potential |2J. 

Here I show that there is a generalized virial theorem for finite scattering 
lengths. I also consider a somewhat more general confinement potential, 

V(r) = r^/(r) 

satisfying (3 > —2, (3 ^ 0, and /3/(r) > 0, where /(f) is any smooth function 
of the unit direction vector f. For a harmonic trap (3 = 1. 

Such a generalized virial theorem is 

E -P±l Ev = -^ (2, 



sir am 



where m is the fermion mass, X = lirm : _ +0 o /c 4 pk<r, and pko- is the momentum 
distribution at momentum Kk and spin state a. The amplitude of pko- is defined 
by / ^yg-pko- = N a , the number of spin-a fermions. Note that iVj and N± are 
arbitrary and may be different. 



X equals the quantity QC in Refs. [3,4|. An equivalent definition of X is given 
in Ref. g]: 

X = lim 7r 2 KNk>K, 

K—>oo 

where Nk>K is the expectation of the total number of fermions with momenta 
larger than hK. 

We will give X a short name: total contact, or simply contact. We will give the 
spatial function C(r) introduced in Ref. [3] a name: local contact density. The 
quantity C = Vt^X = Q^ 1 J d 3 rC(r) [3] will be called average contact density 
(over volume fl). For a homogeneous system of volume Q, C(r) equals C 

To prove (J2j), we first consider an energy eigenstate at scattering length a, 
with energy = E- m t eina \+Ey, where -External is the internal energy expectation 
value. We then modify this state infinitesimally, in two consecutive steps. 

In the first step, we adiabatically change the scattering length to a' = (l + e)o, 
where e is an infinitesimal number. The energy changes to E' = -Ei nternal + E' v . 
Using the adiabatic sweep theorem of Ref. [1] , we find 

h 2 X 

E'-E = -^-e + 0(e 2 ). 



In the second step, we do a geometric compression of the system's wave func- 
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tion, from 0'(ri, ■ ■ ■ , rjy) to 

0"( ri , ■ • ■ , r N ) = (1 + e) 3Ar / 2 0'((l + e)r a , ■ • • , (1 + c)r ff ). 

Using the short-range boundary condition for the wave function (0 oc 1/s — 
1/scatt. length when the distance s between two fermions in different spin 
states is small), we find that 0" corresponds to a state with scattering length 
a" = a' /{I + e) = a. Using the energy theorem [3], we get -Eternal — (1 + 
e) 2 i4 ternal , and E% = (1 + e)^E' v . So 

E" -E' = 2eE i ' nternal - /3eE' v + 0(e 2 ) = 2e£ internal - (3eE v + 0(e 2 ). 



Because the state <\>" has the same scattering length as the initial state 0, and 
because the difference between the two states is of the order e, the variational 
stability of energy levels implies that E"-E = 0(e 2 ), or (E" -E') + (E' -E) = 
0(e 2 ). So 

h 2i? interna i — f3E v = 0. 

Rewriting -External = E - E v , we get $3). 

Obviously, (j2J) is also valid for any statistical ensemble of energy levels, with a 
statistical weight decaying sufficiently fast at large energy, such that X equals 
the statistical average of the values of T's for the individual energy levels [3]. 
Thus 02]) is valid for any finite temperature states in the canonical or grand 
canonical ensemble, as well as the ground state. 

When kpa — > 0~ (hp fixed), I oc a 2 and Eq. (j2j) reduces to the virial theorem 
for the noninteracting Fermi gas. 

When a = oo and (3 — 2, Eq. (j2J) reduces to Eq. (CD). 

When iVj = iVj_ ^> 1, fc^a — > + and the temperature is zero, the system ap- 
proaches a Bose-Einstein condensate of tightly bound molecules, and (j2J) ap- 
proaches the virial theorem for the Gross-Pitaevskii equation for these bosonic 
molecules [5] with scattering length a m ~ 0.6a [6]. 



2 Pressure Relation 



Suppose that the system is in a cubic box of size L = fi 1 / 3 , and a periodic 
boundary condition is imposed. In the absence of the external potential 

n 2 h2c 

P ~ 3 P * = 12W (3) 
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where pe = E/Q is the average energy density, and C = X/Q is the average 
contact density. Equation ([3]) is valid for any energy eigenstate, or any statis- 
tical ensemble of them, with a statistical weight decaying sufficiently fast at 
large energy, such that X equals the statistical average of the values of X's for 
the individual energy levels [3]. This includes any finite temperature states in 
the canonical or grand canonical ensemble, as well as the ground state. 

In Ref. [I], the pressure relation ([3]) is shown for balanced populations of the 
two spin states: iVf = iVj.. 

Here I point out that ([3]) remains valid even if ^ Aj_. This incorporates 
many interesting possibilities, in particular phase separation between super- 
fluid and normal phase \7\ and, consequently, spontaneous spatial inhomo- 
geneity of the energy density and contact density. 

The general proof of ([3]) is very similar to that of ([2]). Starting with any 
energy eigenstate with scattering length a and energy E, we first increase 
the scattering length adiabatically, from a to a' = (1 + e)a, to get E' = 
E + 4 ^ra e + ^( e2 )- We then do a geometric compression of the wave function, 
after which the scattering length changes back to a, the period of the wave 
function becomes L/(l + e), the energy becomes 

E" = (1 + efE 1 = E' + 2eE' + 0(e 2 ) = E + -^-e + 2eE + 0(e 2 ), 

Anam 

and the quantum state becomes <p". 

If we start from the state <ft, and compress the box by a linear factor (1 + e) 
adiabatically, without changing the scattering length, we will get the same 
final state as 4>" . So the pressure is 

„ , E"-E 2E h 2 X 
P = hm = 1 . 

e-o n - (1 + e)- 3 Vl 3fi \2itamVL 

Although Eq. ([3]) is only exact in the absence of external potential, it is ap- 
proximately valid for each local part of the fermionic cloud in a trap, within 
the local density approximation. In this latter case, pE is replaced by the local 
internal energy density, and C is replaced by the local contact density C(r). 
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